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We discuss aspects of the statistical hadronization model for the production of mesons 
with open and hidden charm in ultra-relativistic nuclear collisions. Emphasis is placed 
on what can be inferred from the dependence of the yield of charmonia on the number of 
participants in the collisions. 

1. Introduction 

Experiments with ultra-relativistic nuclei are performed to produce and study the quark- 
gluon plasma. This new state of matter is predicted to exist at high temperatures and/or 
high baryon densities. Numerical solutions of QCD using lattice techniques imply that 
the critical temperature (at zero baryon density) is about 170 MeV [|TJ. Comprehensive 
surveys of the various experimental approaches on how to produce such matter in nucleus- 
nucleus collisions have been given recently ||-f|. Here we focus on charm production 
and its recent interpretation in terms of a statistical hadronization model |J. Since 
the statistical hadronization model for charm is a consequence of the success of thermal 
model descriptions for hadrons produced in nuclear collisions we first very briefly review 
the present state there. We then present an update of the extension of this model to 
charm production with emphasis on what can be inferred from the dependence of the 
yield of charmonia on the number of participants in the collisions. Remarks on what can 
be expected at collider energies will conclude the manuscript. 

2. Equilibration at the Phase Boundary 

The statistical model used to describe hadron multiplicities is presented in detail in [[?]. 
Like its predecessors presented in ||[|] it is based on the use of a grand canonical ensemble 
to describe the partition function and hence the density of the hadrons under consider- 
ation. Without invoking volume information this model can then be used to describe 
ratios of particle yields. The surprizing result of applying this very simple model to data 
is that ratios of all hadron yields including those involving multi-strange baryons, where 
enhancement factors of more than 1 order of magnitude over what had been measured in 
p+Be collisions are observed, can be described consistently J7j if one assumes a fireball 
with temperature T=168 MeV and baryon chemical potential \x& = 266 MeV. This is 
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demonstrated in Fig. [I]. On the other hand, the observed enhancement, especially for 
multistrange hadrons, cannot currently be understood within any of the hadronic event 
generators |T0|| 
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Figure 1. Comparison of measured particle ratios with predictions of the thermal model. 
For details see text and [j7|- 



The chemical potentials \ib and temperatures T resulting from the such thermal analy- 
ses [0-0 place the systems at chemical freeze-out very close to where we currently believe 
is the phase boundary between plasma and hadrons. This is demonstrated in Fig. 0. The 
freeze-out trajectory (solid curve through the data points) is to guide the eye. 

The closeness of the freeze-out parameters (T,/i B ) to the phase boundary might be the 
clue to the apparent chemical equilibration in the hadronic phase: if the system prior to 
reaching freeze-out was in the partonic (plasma) phase, then hadron production in general 
and strangeness production in particular is determined by larger partonic cross sections 
as well as by hadronization. 

3. Statistical Hadronization of Charm 

Charm quarks are heavy (m c T c ) and thermal production of charm quarks and charmed 
hadrons is not likely in ultra-relativistic nuclear collisions. The situation has been recently 
discussed || with the conclusion that, compared to direct hard production, thermal pro- 
duction of charm quarks can be neglected at SPS energies and is small even at LHC 
energy. However, in the course of these investigations we proposed a new scenario for 
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Figure 2. Phase diagram of hadronic and partonic matter. The hadrochemical freeze- 
out points are determined from thermal model analyses of heavy ion collision data at 
SIS, AGS and SPS energy The hatched region indicates the current expectation for the 
phase boundary based on lattice QCD calculations at /U_b=0. The arrow from chemical 
to thermal freeze-out for the SPS corresponds to isentropic expansion. 



quarkonia and charmed hadron production. We assume || that all cc pairs are produced 
in direct, hard collisions, i.e. in line with the above considerations we neglect thermal 
production. For a description of the hadronization of the c and c quarks, i.e. for the 
determination of the relative yields of charmonia, and charmed mesons and baryons, we 
employ the statistical model, with parameters as determined by the analysis of all other 
hadron yields [0. The picture we have in mind is that all hadrons form within a narrow 
time range at or close to the phase boundary. 

Since the number of directly produced charm quarks deviates from the value deter- 
mined by chemical equilibration, we introduced a charm enhancement factor g c by the 
requirement of charm conservation. This leads to the following relation between hard 
open charm production and thermal production: 



N^ ect = -9cC£ N D i : rm + Kr m ) + 9 2 C (J2 N t rm ) + ■■■■ (!) 

i i 

Consequently, the charm enhancement factor g c can be determined if N^ rect is known. 
Since there are no data yet for open charm production in nuclear collisions we rely, as 
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discussed in ||, on the investigations and model calculations of [11]|, where the informa- 



tion available on charm production in hadron-nucleus collisions is extrapolated to Pb-Pb 
collisions at SPS energy. The number of J /if} mesons is then enhanced relative to the 
thermal model prediction by a factors of g^, i.e. 



Nj/* = 9lNf% m . (2) 

Application of this model to the data for charmonium production is somewhat compli- 
cated since the NA50 collaboration has not provided J/if> yields as a function of centrality 
or number of participants in Pb+Pb collisions at SPS energy. Gosset et al. |12[ have 
provided such an analysis based on preliminary 1995 data and this was the basis of the 
discussion in ref . || . We will, for the following discussion, mainly focus on the dependence 
of J /i[j production on the number of participating nucleons. The shape of this dependence 
can be analyzed from published NA50 data fl3HT5|. The absolute normalization we ob- 



tained by using for N part = 100 that ~Nj/^/N h _ = 1.1-10 -6 JL6|, leading to Nj/^/N part 
1.83 -10~ 6 . These results agree well with a recent analysis |I7[ of NA50 data of 1996 and 
1998 within the framework of ref. |l2j and are displayed in Fig. || We note that the so 
obtained data agree with the analysis reported in |T2| for low centralities but exceed those 
data by about 30 % for N part =350. 

In the following we will discuss these data within the framework of ref. |J, with the 
following modifications. We take into account, in the thermal model description, the full 
set of charmed mesons and baryons. This increases, as mentioned already in |J, the 
(grand-canonical) thermal charm yield by about a factor 2.5, mostly because of the large 
statistical factors of the D* mesons. Since, despite this increase, the number of charm 
pairs is significantly less than 1 per collision, we treat the system, as proposed by Redlich 
18], within canonical thermodynamics, following jl{|. A similar approach was recently 



chosen by |20[ , to investigate limits on open charm production which are imposed within 
the present model by the data on J/if> production. 

Neglecting the (very small) quadratic terms in eq. (1) the statistical hadronization 
model then reads: 

1 T (a ]\Jtherm\ 

jw direct _ j\jtherm 1 1 Kile 1 v oc ) /o\ 

cc ~ 2 Vc oc Io ( gcN thermy W 

Here, j\T*^ erm denotes the total thermal yield per collision of open charm mesons and 
baryons (calculated within the grand-canonical ensemble) and I± and Iq are modified 
Bessel functions. From the properties of the Bessel functions one obtaines lim^oo = 1 
so that the grand-canonical limit is obtained for g c Nf£ erm 3> 1- 

Since quarkonia are hidden charm mesons, eq. (2) is still valid. To set the stage for 
the following discussion we now proceed to evaluate the dependence of J/if> production 
on N part . For all estimates we assume N^ rect cx iV^t ■ Since N*^ erm cx N part we get that: 

^ocj?. (4) 

J- 'part 




Figure 3. Comparison of the dependence of the measured (J/i/j)/N pa rt ratio on the num- 
ber of participating nucleons with the predictions of the of the canonical direct/statistical 
model (solid line). The data are from a reanalysis of [|I3| -|T5|1 and are normalized using 
the (J/ip)/h~ ratio deduced in For details, in particular about the normalization, 



see text. The dashed line is obtained by increasing the direct open charm yield by 50 



With the above assumed N part dependences one can then solve Eq. 3 to determine the 
Npart dependence of g c with the result that 



^ « N par, (5) 
* part 

where a = 2/3 in the grand-canonical limit while a = —2/3 in the canonical limit. 
From the smooth behavior of the Bessel functions it is clear that the coefficient a develops 
smoothly[| and continuously from -2/3 to +2/3 as N*r ec< increases from values C 1 to 
values ^ 1. The numerically obtained dependence of a on ~Nf 5 rect for fixed N^; erm is shown 
in Fig. |j. We conclude from this discussion that, in the range of validity of the statistical 
hadronization model, the centrality dependence of the yield will be directly related to 
overall magnitude of the open charm cross section relative to the thermal open charm 
yield. 

1 Numerical solution indicates a slight N part dependence of a for intermediate values of Nfg rect . 
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Figure 4. Dependence of a on the magnitude of the direct open charm yield. For details 
see text. 



Application of this model to the data, displayed in Fig. now proceeds in the following 
way. We start the calculation at N pari = 350 to avoid that part of the data which is 
determined by fluctuations |H] not contained within this model. Furthermore, predictions 
of the model should only be trusted from about N part > 150 on, where also the ip'/(J/ip) 
ratio is close to the thermal value for Pb+Pb data (see Fig. 2 of ref. 0). In this 
context it should be mentioned that the thermal character of the ip'/(J/ip) ratio was first 
noticed by |22j| . In our approach, ratios for all higher charmonia states including the Xc 
should approach the thermal value from N part > 150 on, implying that for those N part 
values feeding to J/ip should be small. In this picture, there should thus not be different 
"thresholds" for the disappearance of different charmonia. 

We also take the value of direct open charm production from PJll| , appropriately scaled 
to N part = 350, i.e. N^ rect = 0.144. With N^ erm = 0.55 for N part = 350 we get, solving eq. 
(3), a value of g c = 1.43. Using eq. (2) and N*^[ m = 1.84 ■ 10~ 4 we obtain immediately 



that 
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J/i> 



part 



1.08 • 10~ 6 for N 



part 



<J/il> 

350 and a = 



—0.61. This result is shown as the 
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solid line in Fig. || Considering the uncertainties in the normalization of the J/ip and 
open charm yields the agreement is surprizing. We illustrate the importance of accurately 
knowing the open charm yield by the dashed line in Fig. [| Here the open charm cross 
section was increased by 50 % (well within the uncertainty of the pp data), leading to 
g c = 1.78 and a = —0.60. Only a direct measurement of open charm production in 
Pb+Pb collisions can remove this uncertainty of the present analysis. 

Since the present analysis involves absolute densities we would like to point out that, in 
the model of 0, an excluded volume correction is applied. For the parameters chosen in |7j 
this leads to a density reduction by a factor of ft ~ 0.65. This should be taken into account 
when trying to deduce limits on open charm production from a statistical analysis [[23 1 . 
Alternatively, one can determine the effective volume by requiring, as is done here for the 
SPS data analysis, that the total charged particle multiplicity is correctly reproduced by 
the thermal model calculations. 



4. Summary and Outlook 

Hadron production results from central nucleus-nucleus collisions at ultra-relativistic ener- 
gies can be quantitatively understood by assuming that the fireball formed in the collision 
freezes out chemically very near to the phase boundary between quark-gluon plasma and 
hadron gas. Mesons containing heavy (charm) quarks are not thermally produced but 
their yield can be quantitatively explained, at least for central collisions, in the statistical 
hadronization model where all charm quarks are directly produced but all hadrons con- 
taining charm quarks are formed according to thermal phase space. For small values of 
j^dvrect ^ ^ ag wag p j n ted out by [|1^] and recently used by |pCfl , important to use the 



canonical approach for calculating the thermal phase space. 

Predicting absolute yields for the rapidity density of J/ip mesons at collider energies 
requires knowledge of the yield of open charm pairs as well as of the total charged particle 
multiplicity (to fix the volume needed for the calculation of thermal yields). The following 
estimates are for central (N part = 350) Au+Au (RHIC) or Pb+Pb (LHC) collisions. They 
should be considered schematic but illustrative of the qualitatively new features expected 
within the framework presented here. For full RHIC energy N^ rect is expected to be of 
order 1.5 per unit rapidity. For a charged particle rapidity density of dN c h/dr] = 1000 we 
get, using eqs. (2) and (3), that g c = 9.0, a = —0.25, resulting in a J/xjj yield of 0.011 per 
unit rapidity, close to the expected unsuppressed value. For LHC energies, where N^ rect 
approaches 50 per unit of rapidity, the grand-canonical picture should be appropriate, i.e. 
a will be close to 2/3. The actual computation for dN c h/drj = 6300 gives a = 0.62 and 
g c = 41, implying a J/ip yield of about 1.4 per unit rapidity, about a factor of 20 above 
the direct yield! In the grand-canonical limit the expected J/ip yields is proportional to 
(N^ rect ) 2 1 'Nch and can be trivially scaled to other values for charm and charged particle 
production. As already discussed in |J, the statistical hadronization picture should, under 
these circumstances, lead to strongly enhanced quarkonia production even compared to 
values for direct hard scattering production. A similar result has been obtained recently 



24] using a model based on kinetic equations. 
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